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I. INTRODUCTION 

The Poisson equation of Newtonian gravitation, 

V2$(t,x) = 47rGp(t,x) (1) 

is a consequence of the inverse-square force law, which is 
ultimately based on solar-system observations that origi- 
nally led to Kepler's laws of planetary motion. Einstein's 
gravitational field equations have generalized equation 
(1) into a consistent relativistic framework that is in good 
agreement with present solar-system data [1-3]. Never- 
theless, on small laboratory scales, for instance, questions 
remain regarding the validity of the inverse-square law of 
gravitation and hence equation (1); at present, efforts 
continue on resolving such experimental problems [4-7]. 
This paper is about deviations from the inverse-square 
force law on galactic scales in order to resolve the prob- 
lem of the flat rotation curves of spiral galaxies. 

An essential component in the conceptual development 
of general relativity (GR) is the way Lorentz invariance 
is employed to describe what accelerated observers mea- 
sure. Lorentz invariance is a fundamental symmetry and 
refers to measurements of ideal inertial observers that 
move uniformly forever on rectilinear timelike worldlines; 
therefore, an assumption is required to relate these ideal 
inertial observers to actual observers that are all non- 
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inertial (i.e., accelerated). The special theory of relativ- 
ity uses the postulate of locality, namely, the assump- 
tion that an accelerated observer is pointwise inertial. 
The hypothesis of locality is known to be an approxi- 
mation [8, 9]; in fact, its domain of applicability is lim- 
ited to motions with sufficiently low accelerations. The 
locality principle is also an essential ingredient of Ein- 
stein's heuristic principle of equivalence that is the cor- 
nerstone of general relativity. Nonlocal special relativity 
is a generalization of the standard theory that goes be- 
yond the locality postulate and involves a certain average 
over the past worldline of the observer [10]. The prin- 
ciple of equivalence of inertial and gravitational masses 
implies a general connection between inertia and gravi- 
tation; therefore, one would expect that the nonlocality 
of accelerated observers in Minkowski spacetime would 
entail a nonlocal theory of gravitation [11]. However, a 
direct nonlocal generalization of GR has not been possi- 
ble; that is, the highly local nature of Einstein's principle 
of equivalence apparently prevents a straightforward non- 
local generalization of GR. On the other hand, gauge the- 
ories of gravitation are in general less restrictive [12, 13]; 
hence, in principle, a nonlocal generalization of GR can 
be constructed within the gauge approach to gravitation. 
Indeed, in recent papers [14, 15], a nonlocal generaliza- 
tion of Einstein's theory of gravitation has been presented 
on the basis of the teleparallel equivalent of GR [16]. In 
the simplest possibility, nonlocality is introduced via a 
scalar kernel. In this approach to nonlocal gravity, non- 
locality can persist in the Newtonian limit of the theory. 

To arrive at this limit in the linear approximation, it 
has been assumed, in addition, that the scalar kernel 
IC{x, y) is a universal function oi x — y and x is supposed 
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to be in the future of y to maintain causality [14, 15]. 
In this case, the nonlocal gravitational field equations 
reduce to 



(2) 



cf. Eq. (62) of [15]. Here G^jy is the linear Einstein 
tensor, k = SttG/c^ and T^^ is the energy-momentum 
tensor of matter {dj^T^" = 0). In this gravitational 
background test particles and light rays respectively fol- 
low timelike and null geodesies of the metric tensor 
gfj^f = Vt^i' + ^A"'' where 77^1/ is the Minkowski metric 
tensor given by diag(l, —1, —1, —1) and /i^j/ is the linear 
perturbation away from flat spacetime. Greek indices run 
from to 3, while Latin indices run from 1 to 3. 

It is useful to move the nonlocal term to the right side 
of Eq. (2) via the Liouville-Neumann method of succes- 
sive substitutions. Let us introduce iterated kernels /C„ 
given by ICi {x, y) = 1C{x, y) and 



lCn+i{x,y) = I ]C{x,z)lCn{z,y)d'^z . 



(3) 



Inspection of Eq. (3) reveals that in each iterated kernel 
JCn{x,y), with n > 1, causality is preserved so that x is 
in the future of y, but in general /C„(a;,?/) is no longer 
a function of x ^ y. This is therefore the ease for the 
reciprocal kernel TZ{x, y) as well. 



-n{x,y) = ^ICn{x,y) . 



n=l 



Thus Eq. (2) can be written as 



Gf,i,{x) = K 



Tf,^{x)+ / Tl{x,y)Tfj,^{y)d'^y 



(4) 



(5) 



so that the nonlocal theory in this approximation is 
equivalent to GR but with an additional source term. 
In fact, the nonlocal aspect of gravity can appear as 
dark matter given by the integral transform of T^i, by 
the causal reciprocal kernel TZ. In this paper we take 
the view that TZ must be determined from observation; 
for instance, lensing observations of colliding clusters of 
galaxies — such as in the case of the Bullet Cluster [17] — 
could provide clues regarding the nature of the full time- 
dependent reciprocal kernel (cf. Section III). 

In the Newtonian limit (c — > 00), retardation effects 
can be neglected and hence we can assume that each 



iterated kernel in Eq. (3) is proportional to 5{x' 
It then follows from Eq. (3) that 

7^(x,2/) = <5(xO-2;°)9(x-y), 



,^0 



y" 
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where q is the spatial convolution kernel. Using this limit- 
ing form of the kernel in Eq. (5), we find in the Newtonian 
limit the nonlocal Poisson equation [14, 15] 



where the "density of dark matter" pD is given by 

PD [t, x) = y '?(x - y)p(i, y)d^y ■ (8) 

Here g is a universal function that is independent of the 
nature of the source [14, 15]. This simplifying assumption 
is relaxed in Section II, where we discuss the general form 
of the nonlocal kernel in the Newtonian limit. 

This paper is based on the assumption that there is 
no actual dark matter. According to the approximation 
scheme employed in [14, 15], the nonlocal aspect of the 
gravitational interaction acts like dark matter of density 
Pd that is linearly related to the actual matter density 
via the kernel q as in equation (8). Consider, for example, 
the circular motion of stars in the disk of a spiral galaxy 
in connection with the observed fiat rotation curves in 
such galaxies (see, for instance, [18, 19] and references 
therein). At radius r outside the bulge, the Newtonian 
acceleration of gravity for such a star is nearly Vg/r, 
where vq is a constant speed. Poisson's equation then 
implies that the corresponding density of "dark" matter 
must be Wo/(47rGr^). Extending this po to a spherical 
distribution of "dark" matter by assumption, the result 
can be compared with Eq. (8): neglecting the extended 
nature of the galactic bulge and setting p{t, y) = M6{y), 
where M is the effective galactic mass, we find 



g(x 



1 



1 



AttX |x-y|2 ' 



(9) 



where A = GM/vq is of the order of Ikpc. The univer- 
sality of the nonlocal kernel implies that A must be a con- 
stant and hence M oc Vq. The resulting nonlocal modifi- 
cation of Poisson's equation (7)-(9) has been previously 
discussed in connection with the Tohline-Kuhn scheme 
[20-22]. In particular, for a point source p{t, x) = Af(5(x), 
Eqs. (7)-(9) imply that 



, GM GM , /Ixl 



(10) 



VH^AttG [p(t,x) + pD(i,x)] 



(7) 



This coincides with Tohline's original suggestion regard- 
ing a modification of Newton's law of gravitation in order 
to account for the flat rotation curves of spiral galaxies 
[20]. A lucid and enlightening account of the Tohline- 
Kuhn approach is contained in the review paper of Beken- 
stein [22]. 

It is clear from this brief account that, as shown in de- 
tail in [14, 15], the Tohline-Kuhn extension of Newtonian 
gravitation to the realm of galaxies can be naturally em- 
bedded within a nonlocal generalization of GR. However, 
the Tohline-Kuhn scheme disagrees with the empirical 
Tully- Fisher law [23]. The Tully-Fisher relation involves 
a correlation between the luminosity of a spiral galaxy 
and the corresponding asymptotic speed vq. This rela- 
tion, combined with other empirical data regarding mass- 
to- light ratio, roughly favors M oc Wq, instead of M oc Vq 
that follows from the Tohline-Kuhn scheme. Various as- 
pects of this issue have been discussed in [22, 24, 25]. To 
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go beyond the Tohline-Kuhn scheme, we discuss a gener- 
aUzation of Eqs. (7)-(8) m Section II, where the Newto- 
nian hmit of nonlocal gravity is discussed in detail. It is 
hoped that this more general treatment of the Newtonian 
limit could help in the resolution of the discrepancy with 
the TuUy-Fisher law. 



II. NEWTONIAN LIMIT OF NONLOCAL 
GRAVITY 

The linear approximation in nonlocal gravity involves 
a linear perturbation away from Minkowski spacetime. 
Consider a background Minkowski spacetime with global 
inertial coordinates .t" = (rf,x). The gravitational po- 
tentials are given by the tetrad field e^^{x) such that 

e^" = '^u + V^"^ , e'^o.^K-Vo., (11) 



We define the modified field strength ^^^p via 
^tiup = C'^i.p + + ''hp - "^au") 

- Vup (V'.M - i'ap.'^) , (13) 

where ip = rjapil''^^ . Thus both Cf^^p and £pi/p arc an- 
tisymmetric in their first two indices. The coordinate 
components of the metric tensor are given by 

gpiy = ?7a/3ep"ei.'^ = ?;^^ + i'^^ + V'^/p ■ (14) 

The nonlocal gravitational field equations in this linear 
approximation scheme then reduce to [14, 15] 

Gp.ix) + v"^ J ^^^^-y^ €,Uy)d'y = «;T,^(x) , (15) 

where G^i, is the linearized Einstein tensor in terms of 
gp.v,T^i, is the energy- momentum tensor of mass-energy 
in Minkowski spacetime and ^^T"^ = 0. Equation (15) 
corresponds to Eq. (60) of [15]. If K{x,y) = IC{x,y), 
then Eq. (15) reduces to Eq. (2), since G'^" d^'t'""' ■ 
In this case, T^i, must be symmetric and hence the an- 
tisymmetric part of ■0^1/ docs not participate in gravita- 
tional dynamics and can thus be neglected. In general, 
however, the energy-momentum tensor is not symmetric; 
therefore, Eq. (15) contains sixteen field equations for the 
sixteen components of 

The kernel K{x, y) could in general depend upon the 
structure of the source in a manner that is consistent 
with the linear approximation scheme. This situation is 
considered in this section for the nonlocal modification of 
Newtonian gravity. However, in previous work [14, 15], 
this possibility was neglected for the sake of simplicity 



and it was assumed instead that K(x, y) is some universal 
function of x — y. 

To approach the Newtonian limit of the nonlocal the- 
ory, we tentatively assume that the dynamics in Eq. (15) 
is dominated by the Newtonian potential $ such that 

l/'OO = "011 = "022 = i'33 = \^ (16) 

and t/j = — 2$/c^, while the other components of ip^^ 
may be neglected in the Newtonian limit. This assump- 
tion corresponds to the circumstance that as in GR, one 
expects that the main effects would be associated with a 
diagonal spacetime metric of the form g^^ = 77^;^ -I- 
where /i^^ = 2c^^$diag(l, 1, 1, 1). In this case, we find 
from Eq. (13) that 



(17) 



moreover, as in GR, c^Gqq — 2V^$. Hence, with Too 



where ippv is proportional to G/c^. The indices are raised pc and 
and lowered by means of the Minkowski metric tensor 
r]ap ■ The gravitational field strength is then given by 



K(x,y)^S{x'' ~y°)k{^,y), 



(18) 



where retardation effects have been neglected, the New- 
tonian limit of Eq. (15) is of the form 



V'$(x) 



E 



aA:(x,y) 9$(y^ 3 



dx'^ 



dy' 



d-'y ^ inGpix) . (19) 



This is a more general form of equations (7) and (8); 
furthermore, for simplicity we have suppressed any tem- 
poral dependence. Equation (19) reduces to Eqs. (7)-(8) 
if A:(x,y) = fe'(x — y), for which the reciprocal kernel is 
q{x — y); however, as pointed out in [14, 15], the kernel 
could in general depend upon the Weitzenbock invariants 
at X and y. For the case under consideration, these are 



— "in"^ 
pvp. — 



G^^G^p" = 9 



9$ 9$ 

dx^dx^ ' 

dx^dx^ ' 
2 



9$ 



my 



(20) 
(21) 

(22) 



It follows that in the limiting case (c — 00) under con- 
sideration, the Weitzenbock invariants all reduce to the 
square of |V$|, which is the magnitude of the Newto- 
nian gravitational acceleration. Hence, one may express 
the kernel as 



fc(x,y) = fc'(x-y) + fc"(x-y; '^"^l 



(23) 



so that fc(x, y) depends on the structure of the source, 
but is otherwise consistent with the linear approximation 
scheme. 

Equations (19) and (23) imply that 



V2$ + J fc'(x,y)V2$d32/ = 47rG(p + p$), 



(24) 
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where p$ is defined by 

p$(x) = --— — > —-—d-^y. (25) 
^ 47rG ^ 7 9a;* 

We recall that (7(x— y) is reciprocal to fc'(x— y); therefore, 

V^* = AnG{p+p^+ J <z(x-y)[p(y)+p$(y)]d3y} . (26) 

In the absence of p$, Eq. (26) is equivalent to equations 
(7) and (8). However, Eqs. (24)-(26) contain a more 
general treatment of the Newtonian limit of the nonlo- 
cal theory. Such a treatment is necessary in order to 
help resolve observational problems associated with the 
empirical TuUy- Fisher relation [23]. 

Equation (26) is a nonlinear intcgro-differential rela- 
tion for the Newtonian potential $. It is clear from Eq. 
(23) that scaling $ by a constant factor leaves the kernel 
invariant. Thus $ given by Eq. (26), despite the nonlin- 
earity of this equation, will be linear in the gravitational 
constant G, as would be expected on physical grounds. 
Moreover, as in Newton's theory, the potential $ can be 
determined from the modified Poisson equation only up 
to an additive constant. Solutions of equations (19) and 
(23), or equivalently Eq. (26), are not known at present; 
therefore, in the following sections we resort to the dis- 
cussion of the solutions of the linear part of the modified 
Poisson equation. 

III. ORIGIN OF KERNEL q 

We return to the study of equations (7) and (8). The 
main feature of these equations is the existence of a linear 
relation between the potential $ and matter density p; 
that is, 

^{t,,^) = G j x{^,y)p{t.y)d^v- (27) 

The Green function x can in this case be simply obtained 
from Eq. (10), namely, x is a function of |x — y| and is 
given by 



x(x,y) = - 1 



1 



(28) 



It can be easily verified that x is a solution of 

V2x(x, y) = 47r [5(x - y) + q(x - y)] . (29) 

One can develop potential theory (sec, for instance, [26]) 
for nonlocal gravity on the basis of equations (27)-(28). 
Moreover, the force of gravity per unit test mass is given 
by 



-V$ -G 



X- y 

|x-y|3 



1 x-y 

A|x-v|2 



pit,y)d^y. (30) 



The integral form of equation (7) — as well as its gener- 
alization in equation (26) — can be obtained using Green's 
theorem; this is the subject of Appendix A. 

In this paper, we take the tentative view that q(r) 
must ultimately be determined via observation. That 
is, this nonlocal "Newtonian" aspect of gravity, just as 
the local Newtonian inverse-square force law, is a fea- 
ture of the gravitational interaction deducible from ex- 
perience. Thus there is no fundamental basis at present 
for the determination of the specific form of the non- 
local kernel other than the concordance of equation (7) 
with observational data. As pointed out in [15], the con- 
volution theorem for Fourier integrals may be employed 
to determine q using Eq. (8) once p and pu are com- 
pletely known. However, this expectation is unrealistic 
at present. Since Newton's time, various modifications 
of the inverse-square force law have been contemplated 
[27]; similarly, we can investigate how the potential (10) 
would change if the kernel (9) is modified. 

Let us first consider a kernel of the form 

g(r) ^ 7-T 2 ^ «2 ' (31) 



AttX r^ 



f2 ' 



where r ~ |r| and £q is a constant length parame- 
ter such that for £o 0, Eq. (31) is, unlike Eq. (9), 
singularity- free. Integrating the corresponding Eq. (29), 
it is straightforward to show that the analog of Eq. (10) 
is in this case 



GM GM 



X 



In 



(r2+ £2)1/2 



A 



1^0 , -1 

— tan 



r V4 

(32) 

We note that the 
ln(£o/A) for r — 



which reduces to Eq. (10) for £q = 0. 
term in square brackets goes to 1 
0; therefore, the logarithmic singularity in Eq. (10) is 
avoided by the introduction of £o ^ 0. Next, let 

1 1 



q{r) = 



AttX r^ 



-r/Lo 



(33) 



where Lq is a constant length that renders the integral 
of Eq. (33) finite over all space. This is necessary to en- 
sure that the total mass of "dark" matter is finite; as 
pointed out in [15], the total "dark matter mass" is infi- 
nite if equation (9) is taken to be valid for |x — y| — > oo. 
But empirical data are not available beyond galaxy clus- 
ters and it is rather likely that Eq. (9) must be modified 
for sufficiently large |x — y|. As before, it is possible to 
integrate equation (29) in this case and the result is 



GM GM 



X 



1 



Lo 
r 



r/L„ 



- 1 



Ei 



(34) 



which reduces to Eq. (10) for Lq = oo. Here G = 
0.577 ... is the Euler constant and we use for the ex- 
ponential integral function Ei (a;) the expression 



Ei(-a;) = C + lnx 



E 



(-l)"x" 



n ■ nl 



x>0 



(35) 
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and the asymptotic expansion 



Ei(-x) 



n=l 



oo ; 



(36) 



see the first formula in (8.214) and formula (8.125) on 
page 927 of Ref. [28]. It is important to note that for 
r — oo, $ has a constant value in this case. 

It is clear from these considerations that variations in 
the simple form of the kernel (9) can lead to complicated 
expressions for the gravitational potential. It is there- 
fore interesting to consider possible unique characteriza- 
tions of this kernel. In this connection, we recall that in 
Newton's theory, the exterior gravitational potential of a 
point mass is proportional to 1/r and satisfies Laplace's 
equation. In fact, the fundamental harmonic solution of 
Laplace's equation in n-dimensional Euclidean space is 
l/r"~^ for n > 2 and Inr for n = 2. In the n = 4 case, 
this result has a natural analog in Minkowski spacetime 
with inertial coordinates {ct,x,y, z), namely, OW = 0, 
where 



w- 



-c''it-tof + {x^xof + {y-yof + iz-zof ■ (37) 



Here □ is the d'Alembertian operator defined by □ = 
—rj'^^dadp . For an interesting discussion of such solutions 
and their singularities, see chapter IX of Synge [29]. 
Let us note, for instance, that Eq. (9) satisfies 



(38) 



That is, up to a constant factor, q(r) is a time- 
independent solution of the scmilinear wave equation [30] 



□ </5 ^ V 



(39) 



It is demonstrated in Appendix B that there is a one- 
parameter family of nonzero spherically symmetric solu- 
tions of Eq. (38) that vanishes together with all of their 
derivatives as ?' — > oo. These solutions, as discussed in 
detail in Appendix B, behave as 



1 

4^ 



1 ^ Co 

^2 'J^^~\~(7 



(40) 



for r — >■ 00 , where a = (vT7— 3)/2 and Co is an arbitrary 
constant parameter. With a suitable choice of Co, these 
latter solutions would also be consistent with galactic 
data. 

Equations (38) and (39) make it possible to contem- 
plate appropriate generalizations of Eq. (9). For instance, 
the invariance of Eq. (38) under spatial translations in- 
dicates that 



1 



1 



AttX {x - xo)2 + {y - yo)'^ + - 



(41) 



is also a solution of Eq. (38) that reduces to q(r),r = 
{x,y,z), for To = 0. Furthermore, it follows in a sim- 
ilar way from the scalar field equation (39) that time- 
dependent kernels can be constructed via Lorentz trans- 
formations. Consider, for instance, a pure boost in the x 



direction with speed v; then, 

1 1 



47rA j^{x + ut)2 4- y2 _|_ 



(42) 



is a solution of Eq. (39) that reduces to q{r) for v = 0. 
Here 7 is the Lorentz factor corresponding to speed v. 
This means that one could construct reciprocal kernels 
involving two events (ci,r) and (ct',r') using functions 
of the form 



1 



1 



AttX 72[(x - x') + v{t - t')]2 + {y- y'Y + [z - z'Y 

(43) 

together with an appropriate causal ordering of the 
events. 

In the rest of the paper, we simply employ kernel (9). 



IV. EXTENDED SPHERICAL SOURCE 

According to the inverse-square force law, a homoge- 
neous spherical distribution of matter attracts an exter- 
nal particle as if the mass of the sphere were concentrated 
at its center. However, this important result of Newto- 
nian gravitation would no longer hold in general with 
a modified force law. To illustrate this point, equations 
(27) and (28) can be used to evaluate the ea;ienor gravita- 
tional potential for any source distribution. Specifically, 
let us consider a spherically symmetric mass distribution 
of radius Rq such that 



M = 4tt I p{r)r dr . 
'0 



(44) 



At a spacetime position (t, X) exterior to the static 
source, i? = jXj > i?o. 



<i>(t,X) 



CM 27rC 

+ 



R 



A 



F{R,r)p{r)r'^dr . (45) 



Here, the first term is due to the fact that in Newto- 
nian gravitation the exterior potential of any spherically 
symmetric distribution can be replaced at its center by a 
point source, whose mass is equal to the total mass of the 
spherical distribution. Moreover, in Eq. (45), F{R,r) is 
given by 



F{R,r)= ( 

JQ 



In 



(i?2 + y,2 _ 2i?r cos 61)1/2 



A 



smf 



(46) 

since the symmetry of the configuration makes it possible 
to choose the z axis to be along the position vector X. It 
is straightforward to show that 



FiR,r)^2\n[j)+fi^), 



(47) 



where for e = r/ R < 1, 



fie) = -1 + _ [(1 + ef ln(l + e) - (l - ef ln(l - e)] 



(48) 
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Thus Eq. (45) can be written as 
GM 



$(t,X) 



R 
27rG 



GAf [R 



A 



f{-)p{rydT, (49) 



which is the sum of the contribution of a point mass M 
as in Eq. (10) and an extra term due to the extension of 
the source. One can show that 



f{e) = -e' + -e^ + 0{e') (50) 
using the following relation that is valid for |a;| < 1, 

ln(l + x) = 5](-l)"+i-. (51) 



To get an explicit result, let us assume, for the sake of 
simplicity, that p{r) = is a. constant. Then, the inte- 
gral in Eq. (49) can be evaluated analytically using the 
formulas (2.729) on page 205 of Ref. [28]. In any case, 
the dominant terms can also be calculated directly from 
Eq. (50) and the end result is 



$ = - 



GM GM 



R 

GM 
lOA 



A 



In 



Rq 



^+14 



i?0 



(52) 



The force of gravity per unit test mass is conservative 
and is given by — V$; for Eq. (52), this points in the 
direction of the source and has a magnitude 



GM GM 



i?2 



\R 



1-if^ 
5 V i? 



1 

35 



R 



(53) 

Here, the quantity in brackets is close to unity, since R > 
Rq; therefore, we may conclude that the extended form 
of a nearly homogeneous spherical source does not signifi- 
cantly alter the main physical results of the Tohline-Kuhn 
scheme. 



term in Eq. (54) is expected to be a very small pertur- 
bation of the Newtonian potential. The following pre- 
liminary considerations are based on the fact that in the 
nonlocal generalization of GR under consideration here, 
light rays and test particles move along null and timelike 
geodesies, respectively. 



A. Time Delay 

In Newtonian gravity the potential vanishes at infin- 
ity by convention. However, the logarithmic term in Eq. 
(54) is assumed to vanish at the radial distance A'. While 
this is of no consequence in Newtonian gravity, it mat- 
ters here as the spacetime interval depends on $. Con- 
sider, for instance, the gravitational time delay A be- 
tween events Pi : (rfi,ri) and P2 : (0^2, r2) when a light 
signal travels from Pi to P2- Let L = |r2 — ri| and 
£ : L he the distance along a straight line from Pi 
to P2; then, A = t2 — {ti + L/c) is given by [31] 



A = — ^ 



Mi. 



(55) 



Pi 



It follows from a detailed calculation that 
n • r2 



2GAf, r2 
— 5 — m — 

ri 

2GM 



n • ri 



c3A 



{(n- r2)ln^ - (n • n) In ^ L 



Here 



-A 



n 



tan 



n • r2 

A 



tan 



n • ri 

A 



^(r2-ri) 



A = -|ri X r2| 



(56) 



(57) 



The net result is a sum of the Shapiro time delay to- 
gether with the contribution of the logarithmic term in 
the potential. The relative magnitude of these delays is 
expected to be similar as in Eq. (54). 

In the rest of this section, only the derivative of the 
potential is involved; therefore. A' drops out of our cal- 
culations. For instance, in the gravitational shift of the 
frequency of light only the difference in the potential $ 
at two spatially separated events is significant. 



V. SOLAR-SYSTEM EFFECTS 

It is interesting to search for evidence of nonlocal grav- 
ity in the solar system. As a first step in this endeavor, 
let us consider gravitational physics in the solar system 
using a Tohline-Kuhn gravitational potential of the form 



GM GM ^ r 



(54) 



where the gravitational source is at the origin of coordi- 
nates. Here A' is assumed to be a galactic-scale length. 
For A = 10 kpc, 2 A.U./A « 10~^, hence the logarithmic 



B. Deflection of Light 

The net deflection angle D of a light ray due to a point 
mass M with potential $ is given by twice the Newtonian 
expectation in the first post-Newtonian approximation, 
which for Eq. (54) works out to be 



D 



AGM 2TrGM 



c2A 



(58) 



where ( is the distance of the closest approach. The 
bending angle is thus slightly larger than the Einstein 
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angle by a constant. However, the extra deflection is not 
expected to remain constant for an extended source (cf. 
Section IV). 

The effect of the logarithmic potential is ~ 10""'^^ of 
the Einstein angle for the bending of light by the Sun. 



C. Perihelion Precession 

The gravitational force due to potential (54) is radial 
and conservative. Therefore, the perturbing influence of 
the logarithmic term in Eq. (54) on Kcplerian orbits is 
such that the orbit remains planar and the orbital angu- 
lar momentum is unchanged. Let (r, 0) be polar coordi- 
nates in the orbital plane and consider an unperturbed 
Keplerian ellipse given by 



where 



a(l-e') 
1 + e cos 6 



(59) 



where a is the scmimajor axis of the ellipse, e is its eccen- 
tricity, (j) ^ (j) — g and g is the argument of the pericen- 
ter. Under the influence of the radial perturbing accelera- 
tion — {GM/X)r~^, the orbital elements of the osculating 
ellipse vary in accordance with the Lagrange planetary 
equations [32]. In this case, we find 



da 2uia?e 
dt^^X{l- e2)3/2 
dg uja 

dt "Ae(l -e2)i/ 



(1 



e cos sm ( 



ecosffl cos< 



(60) 
(61) 



where a; is the Keplerian frequency of the osculating el- 
lipse (w^ = GM/a?). Moreover, [GMa{l - e2)]i/2 jg the 
magnitude of the specific orbital angular momentum and 
remains constant. Let us note here that only positive 
square roots are considered in this paper. 

The elements of the osculating ellipse change slowly 
according to Eqs. (60) and (61); therefore, it is natural 
to average the right-hand sides of these equations over 
the fast orbital motion with period T = I-k juj. That is, 
we define the average of a quantity Q to be 



so that 



< CJ > = 



(1 - 6^)3/2 

27r 



2-K 



Qdt, 



Q< 



(l + ecos( 



(62) 



(63) 



since r'^d(/)/dt = [GMa{l — e^)]^/^ for the unperturbed 
orbit. It follows that < da/dt > = 0, so that the semi- 
major axis remains unchanged on the average. This is 
also the case for the orbital eccentricity due to the con- 
stancy of the angular momentum. Thus the ellipse keeps 
its shape on the average but precesses, since 



dq 



2A 



(64) 



P(e) = 4 f(l - e2)V2 - (1 - e^) 



(65) 



decreases from unity at e 
have used the fact that 



to zero at e = 1. Here we 



da 



27r 



1 -I- ecosa (1 — e^y/'^ 



(66) 



One can obtain the same result for the pericenter pre- 
cession frequency from the study of the variation of the 
Runge-Lenz vector [33]. 

For the solar system, the resulting perihelion preces- 
sion [34] is retrograde and for A = 10 kpc, it is about 
10"'^ of Einstein's value for Mercury and about 2 x 10~^ 
for Earth, as there is more "dark matter" to influence 
the outer orbits. The general relativistic contribution to 
the perihelion precession of Mercury is known at present 
at the level of about one part in a thousand; therefore, 
the possible contribution of the logarithmic potential is 
hidden within the present measurement error. Future 
improvements in such measurements may make it possi- 
ble to detect the influence of nonlocal gravity in the solar 
system. Nevertheless nonlocal effects appear at present 
to be too small to be detectable. For instance, the con- 
tribution of the logarithmic term to the Pioneer anomaly 
is ~ 10~^ of the anomalous acceleration of the Pioneer 
spacecraft. 

It appears that other anomalies in the solar system — 
such as the flyby anomaly, the possible secular increase 
in the Astronomical Unit and the increase in the eccen- 
tricity of Moon's orbit (see Ref. [35]) — are not directly af- 
fected by the conservative perturbing force under consid- 
eration here. We should also mention that solar-system 
deviations from GR can in principle be used to place 
lower bounds on the constant lengthscales that appear 
in the logarithmic term in Eq. (54). 



VI. DISCUSSION 

Starting from first principles, arguments have been ad- 
vanced for a nonlocal generalization of Einstein's the- 
ory of gravitation [10, 11, 14, 15]. In such a theory, the 
gravitational field is local, but satisfies nonlocal integro- 
differential field equations. These are obtained from the 
local field equations via a nonlocal "constitutive" ansatz, 
as described in Appendix C within the general context 
of gauge theories of gravitation that are less restrictive 
than GR and thus make it possible to implement this 
procedure. 

The Newtonian limit of the simplest nonlocal GR the- 
ory involving a scalar constitutive kernel [14, 15] is stud- 
ied in this paper. It is shown that the theory reduces to 
a nonlinear and nonlocal modification of Poisson's equa- 
tion of Newtonian gravity. The exploration of the non- 
linear aspects of this equation is beyond the scope of the 
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present work; therefore, we ignore the nonhnear part of 
this equation and concentrate on the simpler case of the 
hnear Poisson equation. This turns out to be equiva- 
lent to the Tohline-Kuhn scheme of modified Newtonian 
gravity as an alternative to dark matter [20-22, 24, 25]. 
Indeed, on galactic scales, the nonlocal deviation of the 
gravitational interaction from the inverse-square force 
law could be responsible for observational data that have 
been attributed to the presence of dark matter. As a pre- 
liminary step, we study some of the implications of the 
linear nonlocal theory for observations within the solar 
system. 



Appendix A: Integral form of Poisson's equation 

Consider a source density p(x) with compact support 
in a spatial volume V bounded by the surface S. Using 
Eq. (7) and 



V^i — - — r = -47r(5(x - x') 
|x-x'| 

in Green's theorem, we find 

$(x) = -GI + 5, xeV", 



while 



Here 



GI = S, TiiV . 



p(x') + Pd(x')^3^^ 



and S is the surface integral 
1 



S 



47r 



n dn' 'dn'KTZ 



(Al) 

(A2) 
(A3) 

(A4) 
dS , (A5) 



where 7^ = |x - x'|, d<^/dn' := (Vx'^") • n' and n' is the 
unit vector normal to the boundary surface S. 

In Eq. (A3), where x ^ F, and assuming that q and p 
are continuous functions, the order of integration in T — 
when Eq. (8) is taken into account — may be interchanged 
such that 



1= [ /v(x,x')p(x')d^x' 
Jv 



where K is given by 



/v(x,x') = 



1 



q{y - yi')(Py 



(A6) 



(A7) 



IV |x-y| 

If (0 + Pd is bounded for small r = |x|, falls off as 
j,-(2+a) ^j^}^ a > for large r and $ — > as r — > oo, 
then 



<I>(x) 



-G 



P(xO +Pd(x') _^3^, 



(AS) 



To satisfy the conditions for Eq. (AS), we note that for 
a smooth source density p(x) with compact support, pn 
must be finite and g(x) must fall off as r~(^+'') with 6 > 
as r — > oo, since it is simple to check that 



4 = 



,.6+2 



(A9) 



Appendix B: Nonuniqueness of q 

We consider the time independent, spherically sym- 
metric solutions of the nonlinear equation 



(Bl) 



in spacetime. Using spherical coordinates, the spheri- 
cally symmetric, time-independent solutions satisfy the 
ordinary differential equation 



r 



(B2) 



In particular, we wish to determine the solutions that are 
C°°-flat as r — > oo. 

The change of variables (see [36, 37]) 



z = r u. 



or, equivalently. 



In) 



1 



z{t) = e '^u(e'^), u{r) = — z(lnr), 



(B3) 



(B4) 



transforms equation (B2) to the autonomous ordinary 
differential equation 

z" - 3z' + 2z^ 2z^, (B5) 

where the prime signifies differentiation with respect to 

T. 

In the phase plane (see Figure 1), the corresponding 
system 



z' = w, w' = 3w — 2z + 2z^ 



(B6) 



has two rest points {z^w) — (0,0) and [z^vS) — (1,0), 
which correspond to the solutions 



u{r) — 0, u(r) — 



(B7) 



of equation (B2); these solutions might also be obtained 
by inspection of this equation. Both of these solutions 
and all of their derivatives vanish in the limit as r — > oo; 
that is, both are C°°-flat as r oo. 

By linearization at the point (0, 0) in the phase plane 
for system (B6), it follows that this point is a hyper- 
bolic source with corresponding spectrum {1,2}. The 
rest point (1,0) is a hyperbolic saddle with spectrum 
{i(3-yT7), 5(3 + \/T7)}. Moreover, the stable manifold 
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-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



FIG. 1: Phase portrait of system (B6), where w is drawn 
versus z. 

of this rest point is tangent to the hne with parametriza- 
tion 

A solution of the system whose orbit hes on this stable 
manifold is asymptotic to a corresponding solution of the 
linearization on the linearized stable manifold (which is 
exactly the line in display (B8)). That is, the asymptotic 
behavior of a solution on the stable manifold is 

where a = ^(vTT — 3) and tq is an arbitrary real num- 
ber. In particular, the asymptotic behavior of the first 
component of a solution on the stable manifold will be 

z(t) - l±e-'"(^-^°), (BIO) 

which corresponds to 

1 e"'^" 

(Bll) 

Because ct > 0, the function u decreases like as 
r — > oo. 

Appendix C: Nonlocal Poincare gauge theory 

In the traditional approach to general relativity (GR), 
the invariance of the spacetime manifold under coordi- 
nate transformations is emphasized. Thus curvilinear 



coordinate systems play a dominant role in GR and ob- 
servers that occupy fixed positions in space in a given 
coordinate system are accelerated in general. The frame 
field constructed from linearly independent vectors that 
are tangent to the curvilinear coordinate lines are called 
natural or coordinate frames. Such a frame, which is not 
necessarily orthonormal, is integrable by definition and 
hence holonomic. The gravitational field equations are 
second-order local partial differential equations (PDEs) 
for the spacetime metric in a given system of coordinates. 

In the gauge approach to gravity, however, one con- 
siders arbitrary frame fields that are in general non- 
integrable, that is, anholonomic. Taking advantage of 
the freedom afforded by the use of anholonomic frames 
and the associated geometric concepts (such as torsion), 
the gravitational field equations take the form of first- 
order local PDEs [12, 13, 38-40]. 

It turns out that one can extend the first-order local 
field equations to nonlocal ones via the introduction of a 
"constitutive" kernel as in the phenomenological electro- 
dynamics of media. In this way, a nonlocal generaliza- 
tion of Einstein's theory of gravitation becomes possible 
by starting with the teleparallel equivalent of GR rather 
than with GR itself [14, 15]. 

In order to be able to construct a nonlocal generaliza- 
tion of Einstein's gravitational theory, the gauge theory 
of translations was recently employed in [14, 15]. Now, 
the gauge theory of translations itself is a somewhat de- 
generate subcase of the gauge theory of the Poincare 
group, the so-called "Poincare gauge theory of gravity" . 
In turn, the question arises whether the nonlocal gen- 
eralization of the translational gauge theory can be ex- 
tended to the Poincare gauge theory. This is, in fact, 
the case. We will follow the method used in [14, 15] for 
the translational gauge theory in the more general case 
of the Poincare gauge theory. We will use the notation 
and conventions of [15]. 

Let the gauge Lagrangian of the underlying Riemann- 
Gartan spacetime depend on coframe e^", torsion T^" 
and curvature Rij"'^ = —Rij^"; that is, £grav = 
'Cgrav(ei",ry",i?y"^). The matter Lagrangian £mat, 
with the matter field(s) VP, is supposed to be minimally 
coupled to the geometry. Then the total Lagrangian 
reads 

Act = /;grav(e.",T,/,i?,,"'3)+£^,t(er, V, AV-) , (CI) 

with the independent field variables e^" (coframe), 
Vi"'^ = —Ti^"" (Lorentz connection) and (matter 
field(s)). With the help of the two excitations 

the two field equations — the results of the variation of 
i^tot with respect to e^" and P^"^ — can be written as 

DjT-C-' a ^ ~ Sq,' , (G3) 

DjW^aP - e-' [aWyiiS] = Tap'' , (G4) 
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where Sq' = SCraa.t/Sei" denotes the canonical energy- 
momentum tensor density of the matter field and r^p^ = 
i^'Cmat/^ri"^ = —Tpa dcnotcs the corresponding canon- 
ical spin (angular momentum) tensor density (note that 
these definitions differ slightly from the ones in Rcf. [16]). 
The energy-momentum tensor density Sq' should not be 
confused with the torsion tensor Tij°'. 

In Eq. (C3), the energy- momentum tensor of the gauge 
fields can be expressed as 



p i — r — T 



]k 



(C5) 



This is the quadratic local Poincare gauge theory. 

For later purposes, it is convenient to express the La- 
grangian in terms of the excitations: 



1 lin 
2 



1_., 1 qu , 



2k 



-A 



l-Tji] 7^. .Q _ 1 T/ i] n. aP 



(C9) 



On the other hand, for Eq. (C4) the spin of the gauge 
fields is given by ia'^''\j\i5]^ which depends only on the 
translational excitation; thus, it is very simple and we 
have already substituted it directly into Eq. (C4). 

This represents the general framework for the Poincare 
gauge theory. We still have to specify the explicit form 
of the gauge Lagrangian. Following the general scheme 
of a Yang-Mills theory, we assume that the Lagrangian is 
local and quadratic in torsion and curvature. We denote 
the three irreducible pieces of the torsion by •^■^^Ty", for 
/ = 1,2,3, and the six irreducible pieces of the curvature 
by (■^^iJy "'^, for K = 1,2, ...,6; details can be found in 
[13, 16]. Then the Lagrangian reads 



loc 

r 



1 

2k' 



1=1 



2C 



K=l 



(C6) 



where k is Einstein's gravitational constant, A is the cos- 
mological constant and ^ is the dimensionless coupling 
constant of "strong gravity" , which is mediated via the 
propagating Lorentz connection. The constants bi and 
ck are dimensionless and should be of order unity. 

We compute the excitations from Eq. (C6) by partial 
differentiation according to Eq. (C2): 



~9 



(C7) 



1=1 



K=l 



lin . . qu . . 



(C8) 



This Lagrangian will also be valid in the nonlocal case. 
We now generalize the local "constitutive relations" (C7) 
and (C8) to nonlocal ones, again as in [14, 15], by using 
an unknown scalar kernel x(x, x') and the world function 
il and its derivatives for transporting tensors from x' to 
x: 



w^k{x) = -^^-^ I uix,x')n'''n"nkk' 



K 



X x{x,x')J2bj^'^T,,^' {x')./^^d^x', (CIO) 
1=1 



(Cll) 



X x(x,x') CK '^''^R,r''''\x')^J^^)d^x', 

K=l 



lin . . qu . . 



(C12) 
(C13) 



The final field equations will not be written down ex- 
plicitly. We find them as follows: We first substitute 
Eqs. (C10)-(C13) into (C9) and (C5) and then into the 
field equations (C3) and (C4); the last step involves the 
substitution of the new Eq. (C5), after Eq. (C9) is in- 
serted, into Eq. (C3). In this way, we have a set of 16 + 
24 intcgro-differential equations in terms of the variables 
e^°',^^°'^ and 
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